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We measure the state dynamics of a tunable anharmonic quantum system, the Josephson phase 
circuit, under the excitation of a frequency-chirped drive. At smaU anharmonicity, the state evolves 
like a wavepacket - a characteristic response in classical oscillators; in this regime we report ex- 
ponentially enhanced lifetimes of highly excited states, held by the drive. At large anharmonicity, 
we observe sharp steps, corresponding to the excitation of discrete energy levels. The continuous 
transition between the two regimes is mapped by measuring the threshold of these two effects. 
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Ever since the laws of quantum mechanics were for- 
mulated, there has been an ongoing effort to explain the 
emergence of classical laws in experimental systems. The 
first explanation by Bohr states that these systems oper- 
ate in the limit of large quantum numbers jl|, in which 
case they may be described by a wavepacket that on the 
average follows the classical equations of motion [2|. In 
addition, coupling to uncontrolled, external degrees of 
freedom (decoherence) , is often related to the emergence 
of classicality [3|. Recent experiments and calculations 
have demonstrated the quantum to classical transition in 
oscillators, via noise saturation at low temperature due 
to zero point fluctuations [J, |5| , and harmonic behavior 
at high temperatures in a cavity-QED system [6|. 

In a classical anharmonic oscillator, such as a pendu- 
lum, the energy expectation can be deterministically in- 
creased to large values if the driving force is frequency- 
chirped and its amplitude is sufficiently large. This phe- 
nomenon is commonly known as autoresonance [7]. The 
physical mechanism behind this effect is adiabatic, non- 
linear phase-locking between the system and the drive, 
yielding a controllable excitation as the system's reso- 
nance frequency follows the drive frequency as a function 
of time. This effect is utilized in a wide variety of systems 
[8|, M J ^^d recently in Josephson-based oscillators [5, JOJ . 
In a quantum anharmonic oscillator, the expected time 
evolution under a similar drive is sequential excitation 
of single ener gy l evels of the system, or "quantum lad- 
der climbing" [11| . In practice, for a given anharmonic- 
ity the drive itself introduces some mixing between the 
energy levels due to power broadening and finite band- 
width, which may wash out ladder climbing and lead to 
a classical behavior in a quantum system [121 |l3|. In 
this letter, we measure the dynamics in these two dis- 
tinct regimes in the same system by varying the drive 
parameters and the system's anharmonicity. 

Our system, the Josephson phase circuit (JPC, see 
Fig. [1^), is a superconducting oscillator with a nonlin- 
ear inductor formed by a Josephson junction. It can be 
described energetically by a double-well potential that 
depends on the phase difference 6 across the junction. 
We tune the potential by means of an external magnetic 
flux bias fl4] to vary the anharmonicity and measure the 
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Figure 1: Operation and measurement of the Josephson phase 
circuit, (a) Schematics of the circuit and the potential energy 
at different operating biases. The potential shape and anhar- 
monicity Pr are set by the current source /;, and the state 
inside the well is controlled by the microwave drive 7^w. (b) 
State measurement. A short pulse 7mcas is applied in the flux 
bias to selectively tunnel excited levels n > k. The average 
phase 5 is then measured with an on-chip SQUID to detect 
tunneling events. To determine the occupation probabilities 
of all the A^ levels, this process is repeated with a series of dif- 
ferent Zmeaa amplitudes |18l |. (c) Time sequence of the chirp 
experiment. The drive amplitude Q, is expressed in units of 
the Rabi frequency, measured on the first transition. 



state. Traditionally, the circuit is operated as a two- 
level system (qubit) [ij, ll5|, or a d- level system (qudit) 
|16| . by localizing the phase 5 in a shallow well where 
there are only a few energy levels. The quantum state of 
these levels is then controlled by applying nearly resonant 
current pulses. Due to the finite coherence time of the 
system, this generally requires the anharmonicity inside 
the well Pr = (/oi ~ /12) //oi (where fij is the transition 
frequency from level i to level j) to be sufficiently large 
\\J\ . In this work, we vary the anharmonicity over a large 
range (0.002 < /3r < 0.03) in order to tune the system 
between the autoresonance and ladder climbing regimes. 
The occupation probabilities are determined by measur- 
ing the amount of tunneling out of the well due to a short 
pulse in the flux bias that adiabatically reduces the po- 
tential barrier (see Fig. [T]d); because of the exponential 
dependence of the tunneling rate on the barrier height, we 



get a high tunnehng contrast between the states |l6l . Il8| . 
Tunnehng events are detected using an on-chip supercon- 
ducting quantum interference device (SQUID) [I3 • The 
experiment is repeated ~ 10^ times to yield the occupa- 
tion probabihty. 

The time sequence of the experiment is sketched in Fig. 
[TJ;. Our system has negative anharmonicity (/12 < /oi)- 
Therefore, we decrease the drive frequency at a constant 
rate a = 2Trdf/dt, starting higher than the first resonance 
(/oi), in accordance with the phase locking condition. 
The chirp is followed by a measurement pulse in the flux 
bias /meas and the escape probability is measured. This 
process is repeated for different measurement amplitudes 
in order to extract the state occupation probabilities P^ 
[18|. We start measuring the dynamics at a large an- 
harmonicity /3r = 0.023. The time evolution is easily 
understood by looking at the dressed energies of the sys- 
tem in the rotating frame [20| (see Fig. [5^). We start 
the chirp in the positive detuned side (/ > /oi), with the 
system initialized at the ground state. As the chirp pro- 
gresses (decreasing detuning), it reaches an avoided-level 
crossing, associated with the first transition, at the fre- 
quency / = /oi . If the chirp rate a is small relative to the 
splitting introduced by the drive, an adiabatic transition 
[21| (Landau-Zener transition) to the 1*** excited level 
occurs. As the chirp continues, the probability of stay- 
ing on the adiabatic branch (ladder climbing) is higher 
than in the previous transition due to the increased en- 
ergy splitting at higher transitions (/ = fi^i+i). Figure 
[2b shows the processed data of P„ vs. time along the 
chirp for the relevant states n. We clearly observe steps 
in the occupation, corresponding to the ladder climbing 
effect. In phase space (see insets of Fig. ^p for Wigner 
distribution calculated from simulation), the phase is de- 
localized during each step, as expected from a Fock-type 
state dV') = I"-))- In between the steps, there is a partial 
localization of the phase due to the interference of two 
such states. The fidelity of each step in the experiment 
(the degree of correspondence with a Fock-type state) 
decreases as the state number n is increased, as a result 
of the chirp time being comparable to the energy decay 
time (Ti) of the first excited state. 

Next, we measure the evolution during a similar chirp 
but at a much smaller anharmonicity - (3^. = 0.002. Low- 
ering the anharmonicity brings about more mixing be- 
tween the levels for a given drive, and may therefore re- 
sult in the simultaneous excitation of many levels. Figure 
[21; shows the measured time evolution under these condi- 
tions. Instead of sharp steps, we notice a broad excitation 
during the chirp, consisting of up to 6 levels. On top of 
that, we observe large amplitude oscillations, as expected 
from autoresonant wavepacket dynamics i22'|. The oscil- 
lations are also seen in phase space simulation (see inset 
of Fig. [2J3) where the phase of the localized distribution 
(crescent shape) oscillates during the chirp. A detailed 
comparison between data and simulation, made without 
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Figure 2: State dynamics during the chirp, (a) Dressed ener- 
gies of the lowest levels in the rotating frame as a function of 
the drive frequency detuning A from the first transition /oi . 
As the chirp progresses (decreasing A), for a sufficiently small 
chirp rate the state remains on the adiabatic branch (solid 
black line), (b) Measured occupation probability (color-scale) 
as a function of time and level number in the ladder climbing 
regime (/3^ = 0.023, a/2-K = 2MHz/ns, Q./2-k = 27 MHz) and 
(c) autoresonance regime [jSr = 0.002, a/27r = lOMHz/ns, 
Q./2-n = f 90 MHz). The detuning scale in (a) and the time 
scale in (b) are bound by the start and the end of the chirp. 
Insets: simulated Wigner distribution at different times along 
the chirp. 



adjustable parameters is shown in |f8| . 

To check the stability of the generated wavepacket at 
small anharmonicity, we fix the amplitude and frequency 
of the drive at the end of the chirp to their final value (il- 
lustrated in Fig. [3^, in the case fihoid = ^, where fihoid is 
the drive amplitude after the chirp). Figure [5)d shows the 
resulting time evolution after the chirp. The phase-locked 
wavepacket is centered around n ~ 18 and is remarkably 
long-lived, despite the short decay time at these highly 
excited levels. We define the locking probability Piocked 
as the probability to be in the phase-locked state, taken 
for this measurement as the integrated probability for 
levels n > 10 [13|, [l8|. The locking probability decays 
non-exponentially with a time constant Tiockcd = 1.4 /zs, 
where Tiockcd is defined as the time it takes for the lock- 



ing probability to decay to half of its initial value. The 
results of this experiment should be contrasted with the 
measurement shown in Fig. [5J:, where r^hoid = 0. In 
this measurement, the energy expectation (proportional 
to the average level number) decays exponentially at 
roughly Ti « 300 ns, consistent with the expected de- 
cay of a wavepacket in a nearly harmonic oscillator [23j . 
In phase space (insets of Fig. [3t) there is a quick (5 ns) 
delocalization into a pattern of circular fringes due to the 
non-negligible anharmonicity. The short lifetime-limited 
dephasing at (n) = 18 smears out this pattern into a ring 
(30 ns) [2J|, shrinking at a constant rate Fi — 1/Ti, as 
expected. When fihoid = ^ (see insets of Fig. |31d), the 
locked population (crescent shape) remains localized, but 
slowly leaks out through the edge to the unlocked state, 
which freely decays as in Fig. [3]:. 

The results are explained within an effective barrier 
model 
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26l |. where, the drive at the end of the chirp 
and the system's anharmonicity form an effective poten- 
tial barrier for the population that is locked by the chirp. 
In this picture, the size of the potential barrier scales as 
the amplitude of the drive. We find from this theory that 
the resulting lifetime of the locked population is given by 
Tlockod oc exp(77r2hoid/27r) 18|, |25|, where the parameter 
T] depends on the system and drive frequencies [iSj. To 
check this model experimentally, we measure the locking 
probability as a function of time after the chirp and of 
drive amplitude. In this measurement (see Fig. [3ji) the 
chirp parameters are fixed, but the drive amplitude at 
the end of the chirp is varied [27|. We find that Tiockcd 
scales exponentially with fihoid, supporting the effective 
barrier picture. The holding lifetime increases by nearly 
two orders of magnitude to more than 10 /is. The factor 
r\ we extract from this data (j] « 26 ns), is in agreement 
with theoretical prediction (r\ k, 30 ns) and simulation 
(77 K, 24 ns) [18|. Note that in this experiment, the chirp 
is used to prepare the initial locked state only. 

The locking probability is directly measured using a 
calibrated measurement pulse. In Fig. [3^, as the drive 
amplitude is increased near the threshold (^thjI'K sa 
30 MHz), the highly excited (phase- locked) levels become 
more populated, as indicated by the increased escape 
probability at smaller measurement amplitudes. To mea- 
sure the locking probability Piockod, we use a measure- 
ment amplitude that causes only the population in the 
upper levels to tunnel out (dashed line). 

Although the state dynamics during the chirp is fun- 
damentally different at large and small anharmonicities, 
it has common features in both regimes. In addition 
to the notable increase of the system's energy at rela- 
tively small drive amplitudes, both autoresonance and 
ladder climbing have a threshold in amplitude for phase- 
locking. While in autoresonance the threshold ampli- 
tude ^th scales as c?l^, in the ladder climbing regime 
Q^th oc a^l"^ ■ The change in scaling indicates a transition 
between the two regimes J12l|. To map the transition, 
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Figure 3: Decay of a wavepacket. (a) Time sequence of the 
decay measurement after the chirp, (b) Measured occupa- 
tion probability (color-scale) as a function of level number 
and time after the chirp shown in Fig. 2c, with Slhoid/27r = 
190 MHz and (c) ilhoid = 0. Insets of (b) and (c) show the 
simulated Wigner plot at different times along the decay, (d) 
Measured locking probability (color-scale) as a function of 
time and amplitude of the drive after the chirp, with con- 
tours corresponding to Piockcd(ihoid, f^hoid) = 0.5, obtained 
from data, theory and simulation, (e) Escape probability 
(color-scale) as a function of measurement amplitude Jmcas 
and drive amplitude SI after a chirp, with Qf/27r = 10 MHz/ns 
and fir ~ 0.0046. To measure the locking probability, an in- 
termediate /nieas is used (dashed line) at the end of the chirp. 



we measure the locking probability as a function of chirp 
rate, drive amplitude and anharmonicity. 

Following Marcus et al. [l^l we plot the results (see 
Fig. 0^) in the dimensionless parameters space, ri/-\/a 
and /3/-\/a, where /? = 2TTf3rfoi is the absolute an- 
harmonicity [281. The measured threshold, defined by 
Piockcdi^/ ^/oi, (3/ y/a) = 0.5, changes scaling (the depen- 
dence of threshold amplitude on chirp rate) at thresholds 
where /3 w fi (blue line). This condition is met when the 
broadening of the first transition (caused by the drive 
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Figure 4: Transition from autoresonance to ladder-climbing, 
(a) Measured locking probability (color-scale) as a function 
of tiie dimensionless chirp parameters Q./yfa and /3/^/q • 
The red and black lines are the theoretical thresholds for 
autoresonance (Sl"^ = 0.82q^/*/3~^^^) and ladder climbing 
{Vl[l = 0.8a^/^) 13:]. The blue line {Q. = /3) marks the sep- 
aration between the quantum and classical regimes [12l |. (b) 
A simulation of the experiment shown in (a) with the same 
parameters, including the effects of decay and measurement 
at different fir ^■ 



amplitude) is comparable to the frequency difference be- 
tween neighboring transitions. This marks the transition 
between the classical and quantum re gim es, where the 
energy levels are mixed or resolved |12l ll3J . For compar- 
ison, the theoretical threshold lines of autoresonance and 
ladder climbing are shown on the same axes in red and 
black respectively. Our data converges to the theoretical 
scaling at the classical limit. At the quantum limit the 
threshold shows slow oscillations as a function of fij ^fa^ 
centered on the theoretical ladder-climbing threshold line 
with superimposed fast oscillations |18|. The slow oscil- 
lations are reproduced by numerical simulation (see Fig. 
|3b) and are the result of multi- level Landau-Zener tunnel- 
ing effects [13| . In the simulation, the amplitude of these 
oscillations decreases at larger fij \foL values, converging 
to the theoretical ladder climbing threshold scaling [13.]. 
In conclusion, the ability to measure the system's dy- 
namics in different regimes relies on the wide-range tun- 
ability of the Josephson phase circuit. This tunability 
opens the possibility of measuring the full state (state 
tomography) of wavepackets in more coherent devices in 
the future. Using chirps, one can then generate and mea- 
sure "cat-states" [3|, within this macroscopic system. In 
the ladder climbing regime, one can use the chirp to gen- 
erate high fidelity \n) states in lifetime-improved devices, 
without the long calibration process that is commonly re- 



quired. This demonstrates the usefulness of chirped drive 
in creating and manipulating quantum states in the tun- 
able Josephson phase circuit, with applications in rapid 
state preparation and measurement. 

This work was supported by ISF grant 1248/10 and 
BSF grant 2008438. 
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This measurement is carried out with a different sample. 
The results are similar in both samples. 
The complete data set is obtained from mea- 
surements at four different anharmonicities (/3r = 
0.026,0.022,0.004,0.002), in order to span a large 
range of P/^/a while keeping the chirp time T — 
2tt {fin — fan) /a shorter than ~ 21/2, avoiding the ef- 
fects of decay. We find that the threshold amplitude is 
affected by the decay for longer chirp times. The finite 
coherence time T2 only weakly affects the threshold in 
our experiment (la ]. 
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Q ! Materials and methods. The Josephson phase circuit [jj used in the experiment has 
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the following design parameters: critical current Iq ^ 1.5 yuA, capacitance C~1.3pF and 
inductance L^940pH. The qubit has a tunable frequency /oi in the 6-9 GHz range J2|. 
During the experiment the device is thermally coupled to the mixing chamber of a dilution 
refrigerator at 30 niK, where thermal excitations of the qubit are negligible. 

We use a custom built arbitrary waveform generator (AWG) having a fast (1 ns time 
resolution), 14-bit digital-to-analog converter to produce both the chirp signal and the niea- 



CN ■ surement pulse. To produce the chirp, we modulate a high-frequency oscillator, having a 

> ; 

^ ! frequency /lo, using an IQ-mixer. The modulation signals, produced by the AWG, are fed 

<N . 

•^ ' into the 1 and Q ports of the IQ-mixer to give yI{tY + Q{tY cos (27r/Lo^ + <P) at its output, 

t~>^ . where = arctan {Q{t)/I(t)). To produce a frequency shift from the high-frequency oscilla- 



tor, we keep the amplitude at the output constant while varying the phase (p linearly in time; 
to produce a chirp, we use an accelerating phase: (p = 27r/ot — at^/2, where /o = /m — /lo. 



K^ ! /in is the initial frequency of the chirp, and a is the chirp rate. 

Cd I To properly measure the locking probability Piocked, it is generally desirable to have the 

maximal possible chirp bandwidth A/ = fm — fan (/fin being the final frequency of the chirp) 
in order to raise the energy expectation of the locked population higher. This leads to a 
better distinguishability between the locked and the unlocked population at the end of the 
chirp and correspondingly to an increased measurement fidelity of the locking probability, 
as illustrated in Fig SI. The AWG's bandwidth limitation results in an error of up to ~10 % 
in Piockcd at large anharmonicity (Fig. la), however it does not affect the threshold position 
in Fig. 4a. We use the maximal bandwidth (600 MHz), varying the modulation frequency 
from 300 MHz to -300 MHz |3|, and setting the oscillator frequency /lq 200 MHz lower than 
the qubit frequency /oi. The additional 100 MHz of bandwidth beyond the qubit frequency 
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Figure SI: Level occupation of a bifurcated state after chirp in simulation, as a function of ampli- 
tude. Chirp parameters: (a) /?,. = 0.023, a/27r = 6MHz/ns and (b) ^r = 0.002, a/2TT = 12MHz/ns. 
A/ = 600 MHz and /oi — /fin = 500 MHz in both. The locked and unlocked populations are dis- 
criminated by a level cutoff ric (dashed lines), which is experimentally realized by a calibrated 
measurement pulse (see dashed line in Fig. 3e). At large anharmonicity relative to the chirp band- 
width (a), the locked and unlocked populations partially overlap, leading to a maximal error of 
~10% in Piockod for the parameters used in the experiments. 

is taken to reduce the sensitivity of the threshold to initial condition ^. 

Data processing. To extract the state occupation probabilities P^, we use the escape 
probabilities vs. measurement amplitude data ("escape curve", -Posc(-^mcas))- We first measure 
the single-level escape curves by preparing the system in an \n) state, and then measuring 
the escape probability as a function of /meas (see Fig. S2a). Once the single- level escape 
curves P^^^{Ij^ea,s) are at hand, we decompose the measured escape curve of an arbitrary 
state into the single-level basis Pe'sc(-^meas) by optimizing the solution P^ to the set of J 
equations -Pcsc(-^mcas) — J^^riPesci-^Leas) ^ whcrc j = 1, .., J. Generating the \n) state becomes 

n 

increasingly difficult at a larger n, due to the short lifetime of excited states. The procedure 
is even more problematic when the anharmonicity /3 is small and longer pulses are required 
to create the target state with reasonable fidelity. In practice, at the small anharmonicity 
regime that is used in the state dynamics measurement (see main paper), where /3/27r = 
18 MHz, it becomes impossible to prepare the system in an \n) state, even for n > 1. Instead, 
we use the first excited state escape curve, shifted by (5Jmeas(^) = -^mcas(O) — /mcas(^) as 
an approximate escape curve. This approximation is supported by WKB calculation (see 
below). To determine the position of the escape curves we use the chirp data itself: for a 
given state, the measured escape curve contains information about the position of the single 



level escape curves. As seen in Fig. S2a, the position of these escape curves (defined as the 
point where the single-level escape curve increases to 0.5 of its maximal value) is determined 
from the positions of the peaks in the derivative dPcsc{Imca.s)/dt. Due to the finite width of 
the single level escape curves, the peak corresponding to a certain level is visible only when 
the level occupation is sufficiently large. To find /mcas('^) for all the relevant levels, we sum 
the derivative over all the times along the chirp, as illustrated in Fig. S2b. The extracted 
-^meas(^) valucs are plotted in Fig. S2c (red circles). 

Simulation. To check the validity of our estimate for the escape curves, we calculate 
them numerically using the WKB approximation of the level dependent tunneling rates [Sfl: 

Tn = fnexp{-2iSn), (1) 

S3 

where S^ = J \pn{S)\ d6 is the action, 6i are the classical turning points defined in Fig. S2d, 

h 

Pn(5) = \/2m (E'n — U{5)) is the momentum, E'n is the energy of the n*'^ level, U{5) is the 

potential energy, f^ is the classical attempt frequency, m = C ($o/27r) is the effective mass 

and $0 = h/2e is a flux quantum. /„ is calculated using the classical oscillation time: 

/n = 1/r , where t = § dt = 2 / ,fy^ . The energies of the system are calculated by 

diagonalizing the system Hamiltonian: 

^ 2e2 rf2 Jo$o ?^ 1 /^^ ^<l>oV (o^ 

^ = "7^^ - ^ ^°^^ + 2L (*- - ^ j ■ (') 

The circuit parameters are found by best-fltting the calculated lowest frequencies /oi and 
/i2 to the measured ones and flxing the number of levels in the well to 50 (the number 

of levels in the well is obtained from extrapolating the experimental points in Fig. S2c to 

n 

-^meas(-Pcsc = 0.5) = 0) [6]. The siuglc-level escape curves are then given by, Pe"c(-^meas) = 
1 — exp(— rn(/mcas)Ai^)) where At is the measurement pulse length. The calculated positions 
of the single-level escape curves are plotted in Fig. S2c (solid blue line). 

We simulate the state dynamics of our A^-level system under a frequency-chirped drive 
by propagating its density matrix p with the time evolution operator U = exp{iH^At). The 
A^-level Hamiltonian is calculated in the rotating frame of the drive, with the rotating wave 
approximation [7| applied: 
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Figure S2: (a) Left axis: Calculated escape curves of single level states, and of the state |^) = 
l/\/5 (|0) + |1) + |2) + |3) + |4)). Right axis: derivative of the escape curve of |^). (b) Left panel: 
Derivative of the escape curve as a function measurement amplitude and time along the chirp 
shown in Fig. 2c. Right panel: Temporal sum of the data shown in the left panel, as a function of 
measurement amplitude, (c) Experimental and calculated positions of the escape curve. The WKB 
curve is calculated from the level dependent tunneling rates, based on the calculated energies using 
the best fitted circuit parameters, (d) Potential energy of the circuit used for WKB calculation, 
with classical turning points. 
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where en,n+i = '^'^ {fn~i,n — /n,n+i) IS the anharmonicity at the n^^ level and A = A(t) = 
27r (/(t) — /oi) is the frequency detuning of the drive and h = 2TTh is Planck's constant. 
The Rabi amplitude Q is taken as a real constant during the chirp, and the detuning is a 
linearly decreasing function starting at +27r-100MHz and ending at -27r-500MHz, as done 
in the experiment. The anharmonicities en,n+i are calculated from the diagonalization of 
the system Hamiltonian (Eq. [2]). The simulation neglects deviations of matrix elements due 
to the drive, beyond the harmonic oscillator approximation. Namely, for m ^ n ± 1, we 
set (n 5 m) = 0, and for m = ra±lweset (n 5 m) = a/ti + 1, y/n . We find for the 
first order matrix elements (m = ra ± 1) a maximal deviation of order ~ 10^^ at the largest 
anharmonicity, and highest states. The second order matrix elements {m = n ±2) have a 
maximal value of order ~ 10~^, relative to the first order term at the same m value. Higher 
order elements are smaller than ~ 10^^. For m = n, the contribution to the energies for the 
range of drive amplitudes used in the experiment is small compared with the rotating frame 
energies. A separate simulation taking into account all the matrix elements, without the 
rotating wave approximation, yields identical results in the simulations shown below (see 
Fig. S3b,d) to within a ~ 10^^ deviation. Decoherence is taken into account using quantum 
operations (SJ for amplitude and phase damping. 

In Fig. S3 we plot the level populations as a function of time during the chirp shown in Fig. 
2, compared with the experimental data. The simulation is calculated with no fit parameters 
and includes the effect of energy and phase damping. The energies at large anharmonicity 
(Fig. S3b) are estimated from spectroscopy data, while for small anharmonicity (Fig. S3d) 
they are extracted from the diagonalization of the system Hamiltonian. The experimental 
data and simulation agree qualitatively in both regimes. At large anharmonicity, the lengths 
of the "steps" are slightly different in the simulation due to the error in determining the bare 
transition frequencies (obtained from high power spectroscopy, where shifts and broadenings 
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Figure S3: (a) Experimental data and (b) simulation of the dynamics experiment at large anhar- 
monicity (/3/27r=158MHz) shown in Fig. 2b. (c) Experimental data and (d) simulation of the 
dynamics experiment at small anharmonicity (/3/27r=:18MHz) shown in Fig. 2b. 

are significant). At small anharmonicity, we see a smearing of the oscillations at higher states. 
This is mainly due to the frequency dependent drive amplitude. In both measurements (and 



simulations), we used a frequency dependent drive which decreases along the chirp as \/n{t) 
(where n{t) is the expected average state number as a function of time) to compensate for 
the increasing drive coupling at higher states which increases the mixing between the levels. 
This, however, does not affect the locking condition which is determined from the drive 
amplitude at the first transition. 

We compute the locking probability in the simulation by defining a cutoff Uc at intermedi- 
ate levels: Piockcd = XI Pnnj where the level population vanishes (see Fig. SI). The results of 

n>nc 

this simulation are shown in Fig. 4b. All the parameters in the simulation (anharmonicities, 
chirp rates, drive amplitudes, and decay times) are those used/measured in the experiment. 
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We find that the simulation reproduces the main features of the experiment: the position of 
the threshold as a function of (5/ y/a and consequently the transition between autoresonance 
and ladder climbing. The simuation (as well as the experiment), displays slow averaged 
features with superimposed fast oscillations. These features represent interference between 
adjacent levels in the driven system. We observe that the fast oscillations strongly depend 
on initial conditions (the distance from resonance), while the averaged slow oscillations are 
characteristic to the transition between the quantum and classical dynamics. 

The coherence time T2 only weakly affects the threshold in our experiment. This claim is 
supported by the fact that the measured threshold follows that of the decoherence free sim- 
ulation despite the chirp time being longer than T2. In this simulation , we see that at small 
anharmonicity, far off-diagonal elements of the density matrix (high-order coherence terms) 
of a high-amplitude phase-locked wavepacket are negligible and the phase space (Wigner) 



representation of the wavepacket is similar to the one calculated for a classical system J9l. Il0|. 
At large anharmonicity, we find that only the first order coherences of the state (pj,j+i terms 
of the density matrix) are non-zero and they are significantly populated for short times 
(compared to the relevant dephasing time), during transitions between neighboring levels. 

Theory of held drive. The locking time Tiockcd = W"^ (where W is the decay rate 
from the locked state of the nonlinear resonance) is calculated by Dykman et al. in the 
framework of quantum activation [ll|. It is shown that in the case of weak damping and at 
low temperatures {k-oT <^ /i/01), the locking time is given by: 

^locked = cexp(r7fi/27r), (4) 

where, 77 ~ 4/a//oi |/fin — /oil /^r? and c is a constant on the order of Ti. This result is valid 
for intermediate drive amplitudes: 



'4|^n-/oi| „ o /r,_ „ I r f I / 4 I/fin - /oi | /^v 

as is the case in our experiment, where these conditions translate to 4MHz<^ 
VL/2TI <3.7GHz. 

In this theory, the dynamics are considered to be classical while the noise is quantum, 
and is associated with zero-point fluctuations. Moreover, the expression for the locking time 
coincides with the classical formula for the escape time [l2|, when the classical temperature 



m 



12| is replaced by an effective temperature, Tgfj = (/i/m /2fcR) cothf/ifm /2/crT]. A more 



intuitive, but equivalent theory for the locking time is given by Dykman et al. [13| where 
the escape time from an effective potential well associated with the phase-locked state is 
calculated. The potential barrier in this case scales as the drive amplitude. 

We find good agreement between the simulation of this experiment at several anhar- 
monicities and the scaling predicted by Eq. |H The theoretical prediction of the factor 77, 
calculated using the experimental parameters (see black dashed line in Fig. 3d) is within 
15 % from that obtained in the simulation with the same parameters. 
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